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i-5 Abstract 

Oh 

+-> The problem of classical particle in linear potential is studied by using the formalism of Hilbert space 

and tomographic probability distribution. The Liouville equation for this problem is solved by finding 
the density matrix satisfying von Newmann-like equation in the form of product of wave functions. 



^h The relation to quantum mechanics is discussed. 
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1 Introduction 

q 

tH- Recently [I] the tomographic probability representation of quantum states was introduced. In the 

probability representation the states are identified with the fair probability distributions connected with 
wave functions by integral transforms [2| [31. The density matrix of mixed quantum states can be also 
j* related to the tomographic probability representation [4] . 

The Wigner functions (5J is determined by the tomographic probability distribution in view of Radon 
transform (6J. In [7] the tomographic probability distribution was introduced to describe the states of 
classical particles. This approach was based on using the Radon transform of the probability distribution 
function f(q,p, t) on the phase space. The tomographic probability was called the state tomogram and its 
properties were discussed in M in the context of studying the classical and quantum chaos. Thus, there 
exists the duality of considering classical and quantum systems. One can formulate quantum mechanics 
by using the mathematical tools of classical statistics like probability distribution functions considered 
as alternative of wave functions. On the other hand, one can introduce the quantumdike formalism 
of Hilbert space, wave functions and density operators considering classical systems. The idea of the 
introducing the Hilbert space formalism for classical system states was suggested in |9|. But in the 
tomographic approach the quantumdike formalism used in classical domain is different from suggested 
by Koopman, though the general idea is related to the approach presented in [9] and other works [To] |11|. 



Some relations of quantumdike and classical descriptions were studied in [12] 13 



Tomographic representation of quantum and classical mechanics was reviewed in [14] . It was shown 
that quantum states with discrete variables like spin also can be described by using probability distribu- 



tion functions called tomograms [15] 16 . For continuous variables we got an idea of describing classical 
states by using formalism of quantum mechanics, but for discrete variables such approach is not devel- 
oped yet. Thus, in our work we concentrate on introducing description of the states of classical systems 



with continuous variables by using formalism of quantum mechanics and following 17 . 

In this work we consider one example of classical systems. It is the problem of a particle moving in 
the linear potential. On this example we consider some connections between the standard description 
and Hilbert space representations of classical mechanics. 

The paper is organized as follows. In next section Sec. 2 we describe our classical system using 
standard approach. Then, in Sec. 3 we introduce density matrix by using the Wigner function and its 
similarity to classical probability distribution function f(q,p,t). In next two sections, Sec. 4 and Sec. 
5, classical wave function for gaussian state is found and also relation between evolution operator of 
Liouville equation and Green's function of wave equation is discussed. The conclusions are presented in 
Sec. 6. 

2 Review of the classical problem 

Let us formulate our problem. We consider classical particle motion in the field with linear potential 
V = q. Here q is a position of the classical particle. We use positive sign for potential. For instance, our 
problem can be considered as the problem of electron's moving between plates of flat condenser. Also 
we take the mass of our particle and its charge m = Q = 1. Thus, we take the Hamiltonian of particle 
moving in linear potential in the form 

So, we get the equations of motions for the position q and the momentum p 



dH . dH 
That yields 



dp ' dq 



•7 + 1 = 0. (2) 

It is easy to solve this equation, so let us write the solution 

P = Po~t, 

t 2 (3) 

q = qo+pot- —■ 

Here po is initial momentum, go is initial position. 

In case position and momentum fluctuate, the state of our system is described by a probability 
distribution function f(q,p,t) on the phase space. This function is nonnegative and it satisfies the 
normalization condition, i. e. 

f(q,p,t)>0, 

(4) 
f(q,p,t)dqdp= 1. 



The function is a solution of the Liouville kinetic equation 

d£ + dHdl _ OHdl = Q 
dt dp dq dq dp 

In case of linear potential the Liouville equation reads 

d JL +p d JL_?L = {) (5) 

dt dq dp 

In next sections we will consider solutions of the above equation by applying the quantum- like formalism 
of the Hilbert space vectors and the corresponding wave functions. 

3 Classical wave function 



So, we want to consider solution of equation (5) by using quantum-like formalism. But before we 
start we should show how we introduce this formalism in classical mechanics. In quantum mechanics, the 
states are described by a wave function or a density matrix. But also we can use other representations. 
For instance there exists the phase-space quasidistribution W(q,p,t) called the Wigner function. It is 



given by the Fourier transform of a density matrix p(x,x',t) introduced in 18 

u u 

r q ~2 



W(q,p,t)= I p(q+^q-^,t)e-^du. (6) 



The inverse Fourier transform yields 



p( x , x ',t) = ±- [w(^-,p,t)e*<*-*)dp. (7) 



2tt 

So we get the relation between this two representations of quantum mechanics. They provide two 
ways to describe quantum states. Doing the transform (6) we associate the quantum state with the 
function on phase space. 

In classical mechanics we are usually (we can say always) work in phase space. It is the standard 
tool. But it is possible to make transform to use the Hilbert space representation by introducing density 
matrix (density operator) given formally by (7). So, let us consider probability distribution f(q,p,t) 
which is an analog of the Wigner function. It means that we replace W(q,p,t) by 2irf(q,p,t) in formula 
(7). Then we can write expression for density matrix of classical motion 

p(x,x',t) = ff(^±^, p ,t) e^-^dp. (8) 



Then it is important to introduce the evolution equation. In our case we get the equation which 
resembles the quantum equation. In fact, doing some transformations with Liouville equation (5), finally 
we get von Neumann-like equation for the density matrix 

dp(x,x',t) 1 / d 2 d 

i 



dt 2 



f d 2 d 2 \ 

\Jh? ~ dx 12 ) p ( x ' x> ' *) + <yX ~ x '^ p ( x ' x '' ') ' (9) 



where we used the substitution rules 



df dp(x,x',t) 

f 
dp 



dt dt 

1 i(x-x')p(x,x't), (10) 



P 8q 1i \dx 2 dx' 2 ) PK ' ,h 

This equation is equation with separable variables x, x' . In view of this one can find the solution of 
the equation factorizing the density matrix in the form 

p(x,x',t) =^(x,t)^* (x',t). (11) 

Then, one can show that the function *$>(x,t), which can be called the wave function of classical particle 
satisfies the Shrodinger-like equation. 

But let us pay our attention on the following circumstance. The formula for the density matrix 
provides the possibility to make the gauge transformation of the wave function 

^(x,t)^^(x,t)e i ^ t \ (12) 

Here <p(t) is a real phase factor depending on time. From this point of view, the classical probability 
density f(q,p,t) can be mapped onto density matrix of the factorized form p(x,x',t) = VP(£, t)^>*(x', t). 
The corresponding wave function is determined up of this gauge factor. But we want the function ^(x, t) 
to satisfy Shrodinger-like equation. So, when we use this condition we get equation for the wave function 

.dV(x,t) ld 2 ^(x,t) T , , 

It is worthy to point on that we lose time-dependent phase-factor when we get Shrodinger-like equation 
from the von Neumann-like equation for density matrix. 

Also we must take into account that not each solution of wave equation provides the real probability 
distribution function. The function f(q,p,t) expressed in terms of solution *&(x, t) of Shrodinger-like 
always satisfies Liouville equation, but this solution can take negative values which is not admissible for 
classical probability distribution. 

4 Gaussian solution 



In this section we will consider gaussian solution of our problem and show some connections between 
these two representations. Firstly we will find solution of the Liouville equation, in another words we will 
get probability density function and using (7) we will find the density matrix. The Liouville equation can 
be reduced, by using the Fourier transform, to the von Neumann-like equation as it was demonstrated. 
Let us suppose that the probability distribution at t = has the Gaussian form 

fo(qo,Po) = f(q,P,0) = -e-rt-P 2 o. (14) 

7T 



We want to find evolution of the probability distribution function. This problem can be solved easily. 
In fact, the equations of motion read 

P = Po-t, 



q = qo+ pot 



t/ 



For the initial momentum and the initial position we find the expressions 

Po = P + t, 



qo = q-pt 



t z 



(15) 



Since the probability density is the integral of motion from the Liouville theorem we get the probability 
distribution function for arbitrary time moment 



f(q,p,t) = -e v 

7T 



if 



(16) 



Also we can find the evolution of statistical parameters of our Gaussian probability distribution like 
variances and covariance of random position and momentum. 
We have the initial dispersion matrix 



C n 



1 



1 
1 



(17) 



The transformation matrix, corresponding to the evolution (16) reads 



A 



1 t 
1 



(18) 



Then we get the dispersion matrix for arbitrary time moment t 



C = AC n A T 



l + t 2 t 



(19) 



Our aim is to find the wave function that is associated with this probability distribution function. 
We put expression for f(q,p,t) (16) in (8), and get the density matrix 



p(x,x ,t) 



1 



7T 



2\ 2 



-(p+*) 2 -(g-p*-V) Jp(x-x') 



dp. 



(20) 



Factorizing the expression for the density matrix p(x,x',t) we get this wave function \l/(x, t) 

1 t 2 /2+2x 2 +2t 2 x+4ixt-2ix 2 



*(x,t) 



cxp 



(T+F) 



lix 2 t+2it 3 x \ 



$/W+v) 



(21) 



There is the problem of the time-dependent phase-factor as pointed out in the end of Sec. 3. That's 
why we have one problem, this wave function (21) doesn't satisfy Shrodinger-like equation. So, we can 
see it 

19 2 .d \, , 1 /4 + 8i 2 + 3i 4 + i 6 \ T , . /n 

w-'ft + *r (a '' )= 8( (i + ^ r^°- ^ 

But we want the wave function of classical particle to satisfy the Shrodinger-like equation. Because 
of the gauge invariance, in general, the wave function has the form 

expl-l t2/2+2x2+2t2x ^~ 2ix2t+2it3x +i(l>(t)) 
* (x, t) = ^-= 7=— -■ (23) 

As we require, this function has to satisfy Shrodinger-like equation (13). Using this condition and 
putting (23) in (13), we find the equation for time dependent phase 4>(t) 

4 + 8t 2 + 3t 4 + t 6 , ■ \ T . . 

8(1 + t 2)2 + * W ) * &> ' ) = - ^ 

Solving the differential equation for the phase (f>(t), we finally get 

t t 3 1 . , 1 / t \ 

<P{t) = -g - 24 ~ -arctan(t) + - [^+1?) + consL v 2 5) 

In our case, for this choice of time-dependent phase, the wave function ^(x, t) satisfies the Shrodinger-like 
equation (13). 

5 Relation between propagator of the Liouville equation and Green 
function of the wave equation 

In previous section we found gaussian solution of our classical problem. We presented two ways to 
describe our classical system. We found probability distribution (16) and corresponding wave function 
(23). Also we can find this wave function using Green function of the Shrodinger-like equation. In this 
section we want get it. 

We are going to find the relation between the propagator of the Liouville equation (5) and the Green 
function of the Shrodinger-like equation (13). 

Let us introduce the following maps: 

P -> /, Po -> /o, /o ->• /, Po ->• P- (26) 

The operators F, F~ l , II, B can be given in matrix form expressed as kernels of the integral transforms. 
These maps are determined by the following expressions 

P(x,y,t)= / F~ 1 (x 1 y 1 q,p,t)f(q,p,t)dqdp 



and 

f{q,P,t)= I F(x,y,q,p,t)p(x,y,t)dxdy, 



here from (6), (7), (8) we find 

' x + y 



F- 1 (x,y,q,p,t) = e^ x -y'>5^ 



*■ -iv(x-v) r I x ' y 



F(x,y, q ,p,t) = ±-e-^-^5( 



2 

f(q,P,t) = / n(q,p,q',p',t)fo(q',p')dq'dp' 
here H(q,p,q' ,p' ,t) is propagator for the Liouville equation, 



p(x, y,t)= / B(x, y,x ,y , t)p (x , y )dx dy , 

the matrix B is expressed in terms of the Shrodinger-like equation B (x, y, x' , y' , t) = G (x, x' , t) G* (y, y' , t) 
Then we get identity 

Bp = F- l tlF Po , (27) 

that provide this expression 

B (x, y, x', y', t) p (x, y) dx'dy = 



1 /• _ f f ,,// 8 

= — F 1 (x,y,q,p,t)dqdp dq'dp'U(q,p,q ,p ,t) / F(x ,y' ,q' ,p ,t)p (x' ,y') dx'dy' . 

We compare both part of this identity and get expression for function B(x, y, x' , y' , t). So, after some 
transformations we find expression for function B(x,y,x' ,y' ,t) 

B (x,y,x',y',t) = L j e **-v) e -*V-lt)n (q,p,q' ,p' ,t) 5 (q - *-±A 5 (q' - *±^\ dqdq'dp'dp. 

(29) 
Propagator for the Liouville equation (5) reads 

R{q,P,q',p',t) = 5 (q - q (q,P,t)) S (p - po(q,p,t)) . (30) 

Then, we get 

/* 1 I \ q ~ q ' ~ ^~) 1 / x + y\ ( x' + y n 

B(x,y,x',y',t) = / — exp I i (x - y) } 5 ( q - —^— ) 8 ( q' - —^ — ) dqdq . (31) 



That provides Green function 



1 . (x-x') 2 -t 2 (x+x') 

G(x, x', t) = -== e 1 " « " . (32) 



In fact, because of gauge invariance Green function has the form 



1 .lx-x') 2 -t 2 (x + x') .,, . 

G(x, x', t) = -L= e*- L ^r A l+ ^ l \ (33) 



Indeed, we remember that we lose the time-dependent phase. We remember our demand that Green 
function like wave function must satisfy the Shrodinger-like equation. So, using last condition we can 
find the time-dependent phase up to a constant. Expression for the time-dependent phase reads 

t 3 
<t>(t) = -— + const. (34) 

Thus we find Green function of the Shrodinger-like equation. It gives us opportunity to get all possible 
wave functions and these functions can describe or determine the classical states. But we should take into 
account the nonnegativity property of probability distribution on phase space. Not each wave function 
provides us with a real probability distribution. Using transform (6) we can find solution of the Liouville 
equation. But possibly this solution will not satisfy conditions (4) (for instance this function can take 
negative values) and that is why it will not be probability distribution function, in another words it will 
not describe real classical state. In general, solution of the Shrodinger equation always provides us with 
solution of the Liouville equation but it does not have to correspond to real classical state. 

6 Conclusions 

To resume we formulate the main results of our work. We have constructed the solution of Liouville 
equation of a classical particle moving in linear potential in the form of Fourier transform of the product 
of two " wave functions" . The wave functions satisfy the Shrodinger-like equations identical to Shrodinger 
equations for quantum particle moving in the linear potential. We demonstrated on this example that in 
classical mechanics one can use formalism of Hilbert spaces and formalism of quantum density operator. 
It is another example in addition to examples considered previously in [9J and [17] which provides the 
possibility to use the same formalism both in classical and quantum mechanics. 
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